Can one measure negative specific heat in the canonical statistical ensemble? 



F. Stanisciai^2^ A. Turchi^, D. FancUi"* P.H. Chavanis^'S, G. De Ninno'^-^ 
1. Dipartimento di Fisica, Universita di Trieste, Italy 
2. Sincrotrone Trieste, S.S. 14 km 163.5, Basovizza (Ts), Italy 
3. Dipartimento di Sistemi e Informatica, Universita di Firenze, via S. Marta 3, 50139 Firenze, Italy 
4- Dipartimento di Energetica "Sergio Stecco", Universita di Firenze, via S. Marta 3, 50139 Firenze, Italy 
5. Universite de Toulouse, UPS, Laboratoire de Physique Theorique (IRSAMC), F-31062 Toulouse, France 
6. CNRS, Laboratoire de Physique Theorique (IRSAMC), F-31062 Toulouse, France 
7. Physics Department, Nova Gorica University, Nova Gorica, Slovenia 
(Dated: March 2, 2010) 

According to thermodynamics, the specific heat of Bohzmannian short-range interacting systems 
is a positive quantity. Less intuitive properties are instead displayed by systems characterized by 
long-range interactions. In that case, the sign of specific heat depends on the considered statistical 
ensemble: negative specific heat can be found in isolated systems, which are studied in the framework 
of the microcanonical ensemble; on the other hand, it is generally recognized that a positive specific 
heat should always be measured in systems in contact with a thermal bath, for which the canonical 
ensemble is the appropriate one. We demonstrate that the latter assumption is not generally true: 
one can in principle measure negative specific heat also in the canonical ensemble if the system 
under scrutiny is non-Boltzmannian and/or out-of-equilibrium. 

PACS numbers: 



In classical thermodynamics, specific lieat is defined as 
the energy required to increase the temperature of a unit 
quantity of a substance by a unit of temperature. Such 
a definition follows common sense: providing energy to a 
physical system should induce heating, that is an increase 
of the system temperature. In fact, this is normally the 
case of physical systems obeying the Boltzmann statistics 
[H and being characterized by "short-range" interactions. 
In a short-range interacting system, the interaction po- 
tential between different components decays as l/r" with 
a > d {r being the distance and d the system dimension) . 
However, less intuitive thermodynamic properties may be 
displayed by "long-range" interacting (LRI) systems, for 
which a < d. In that case, experiments realized on iso- 
lated systems, described by the micocanonical statistical 
ensemble, may give different results from similar exper- 
iments performed on systems in contact with a thermal 
bath; for which the canonical ensemble is the appropriate 
one. According to literature (see, e.g., and references 
therein), one of the most striking features of such ensem- 
ble inequivalence is the presence of a negative specific 
heat in the microcanonical ensemble: increasing the en- 
ergy of an isolated LRI system may lead to a decrease of 
its thermodynamic temperature. On the other hand, it is 
generally recognized that the specific heat in the canon- 
ical ensemble is always positive, even if interactions arc 
long range. Hence, it is a generally shared opinion that 
one can not measure negative specific heat in a system 
in contact with a thermal bath. As we will see, such a 
result relies on the assumption that the thermodynamic 
temperature, defined as Tth ~ {ds/de)~^ (where s and e 
are, respectively, the system entropy and the energy per 
particle), coincides with the kinetic temperature, Tkin, 
defined in terms of the average of the system kinetic en- 
ergy. In fact, this is true only for Boltzmannian systems, 
and generally not for systems following a different statis- 



tics. 

In this Letter, we demonstrate that negative specific 
heat in the statistical canonical ensemble can be mea- 
sured in non-Boltzmannian and/or out-of-equilibrium 
systems. Moreover, the common opinion according to 
which measuring a negative specific heat (in the micro- 
canonical ensemble) is to be considered as a signature of 
ensemble inequivalence should be revised: we will show 
that ensemble equivalence may still hold when negative 
specific heat is measured (in both ensembles). 

As a paradigmatic example for our investigation, we 
consider the so-called Hamiltonian Mean Field (HMF) 
model , which has been widely studied in the past as a 
prototype of LRI system. The HMF model, which shares 
many similarities with gravitational and charged sheet 
models, describes the one-dimensional motion of N ro- 
tators coupled through a mean field cosine interaction. 
The system Hamiltonian reads 

N N 

where 9j represents the orientation of the j-th rotator 
and Pj stands for its conjugated momentum. To monitor 
the system evolution, it is customary to introduce the 
magnetization M , an order parameter defined as 

II V m-l 

M = where = (cos6'i,sin6'i). (2) 

The infinite-range coupling between rotators, provides 
the system dynamics with all typical characteristics of a 
LRI system. In particular, it is well known that after a 
transient regime (the so-called "violent relaxation"), and 
before attaining the final Boltzmann equilibrium, the sys- 
tem may be trapped in a "quasi-stationary state" (QSS), 
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whose hfetime diverges with the system size N 1,13. In 
this regime, the magnetization is lower than predicted by 
the Boltzmann equilibrium and the system displays non 
Gaussian velocity distributions Q. It has been shown 
that QSS's can be related to the stable steady states of 
the Vlasov equation describing the system in the limit 
— )■ cx) 0. The idea was inspired by the seminal work 
of Lynden-Bell , developed in the context of stellar dy- 
namics, and later applied to vortex dynamics (see, e.g., 

i)- 

Lynden-Bell's approach goes as follows. Starting from 
an arbitrary initial condition, the system evolution makes 
the single-particle distribution function stir in phase 
space. However, as a property of the Vlasov equation, 
the hypervolumes of its values 77 (level of phase den- 
sity), arc conserved. Let us introduce the probabil- 
ity density p{9,p,r]) of finding the levels of phase den- 
sity T] in a small neighborhood of the position {0,p) 
in phase space. Using p, Lynden-Bell defined a lo- 
cally averaged ("coarse-grained") distribution function 
f{9,p) = j p{9,p,ri)ridri, and, in order to determine 
the system equilibrium distribution, introduced an en- 
tropy functional, like in ordinary statistical mechanics: 
s = — J p{9,p,r))lii p{9,p,r])dr/d6dp. Following Lynden- 
Bell, we consider the particular case of two-level water- 
bag initial conditions, i.e. ry = /o or = 0. Within this 
approximation, the entropy reduces to 
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where / = fap. As we will show in the following, max- 
imizing this entropy, one obtains a "fermionic" solution, 
/qss(^,p), which represents the most probable distribu- 
tion when the system is trapped in the QSS, and which 
coincides with a particular class of stable stationary so- 
lutions of the Vlasov equation. While QSS's represent 
out-of-equilibrium states of the A^-body dynamics, they 
could be equally interpreted as equilibrium configuration 
of the corresponding continuous description: in this re- 
spect, the conclusions of our analysis will apply to both 
equilibrium and non equilibrium dynamics, provided the 
latter bears distinctive non-Boltzmannian traits. 

Before proceeding further, it is worth emphasizing that 
our study is relevant for a quite large class of systems and 
for many practical physical situations. In fact, experi- 
mental LRI systems often involve a very large number 
of components. When this is the case, QSS's may last 
long enough to prevent the system from reaching the final 
Boltzmann equilibrium during the finite time of an ex- 
perimental session. For all these experiments, QSS's rep- 
resent the effective (measurable) equilibrium states. As 
an example, one can mention charged systems, e.g. plas- 
mas 0] , systems involving wave- (many) particle interac- 
tions, e.g. free-electron lasers and collective atomic 
recoil lasing (CARL) and magnetic dipolar systems 
[1^ . Long lasting QSS's may be also displayed by grav- 
itational systems (e.g., galaxies), to which Lynden-Bell 



theory was first applied. The model we are here consid- 
ering, although very simple, reproduces the most relevant 
features of the above mentioned systems. In particular, it 
can be shown that, under some hypotheses, the Hamilto- 
nian ([!]) can be formally reduced to that governing CARL 
dynamics [Tsj . This important observation will make it 
eventually possible to test our prediction versus direct 
experiments. Further, it is important to stress that for 
systems characterized by small energy dispersions, the 
two-level water-bag initial conditions represent a good 
approximation of a more natural Gaussian initial distri- 
bution [3. 

In the following, we briefly review the standard argu- 
ment leading to the conclusion that the specific heat in 
the canonical ensemble is always positive, and we discuss 
in which respects such an argument does not apply to 
our case. 

Consider a generic (long-range or short-range) sys- 
tem. Given the microcanonical entropy s(e, n), where 
n is the particle density, the canonical rescaled free 
energy, (f){(3,n) = /?/(/?, n), can be calculated as the 
Legendre-Fenchel Transform (LFT) of s(e, n): (j){l3,n) = 
infe [/3e — s(e,n)]. Here /? is a free parameter to be re- 
lated to the system temperature. It is worth reminding 
that the LFT of a generic function is always a concave 
function. If s{e,n) is also concave, the inverse LFT can 
be applied: s(e, n) = inf^ [/3e — (f>{/3,n)]. This is always 
the case of short-range systems, for which the equivalence 
of statistical ensembles holds true. Using the previous re- 
lation one can write 



as 



de 



where 



de 
df 



de^ 



dT 
'd/3 



(4) 



(5) 



is the specific heat (at fixed volume) calculated in the 
microcanonical ensemble and T is either the thermody- 
namic (Tth) or the kinetic (Tkin) system temperature. 

Making use of the LFT, the specific heat can be also 
calculated in the canonical ensemble and expressed in 
terms of the rescaled free energy: 



9^ 



dT 
'dp 



(6) 



Using the thermodynamical definition of tempera- 
ture: Tth = {ds/de)-^ = one gets (c™^)^^^ = 

-l/T^i^ id'^s/de^). The sign of (C''^)^j^ depends on that 
of {d'^s/de^y. for concave entropies (short-range sys- 



tems), i.e. d s/de^ < 0, 



'th 



> 0; if the entropy 



has a "convex intruder", i.e. d s/de > in some energy 
< 0. The latter may be the case of LRI 



range, (c'^ 
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systems in the presence of a first-order phase transition 

0. In the canonical ensemble, {cl^^)^^^ = —^r i^^^- 

As we have already mentioned, is always concave and, 
as a consequence, (c5;™)^j^ is always positive. 

Consider now the definition of kinetic temperature, 
which is the one experimentally accessible 
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tkin 



p^F{e,p)d0dp 



(7) 



(we assume unitary mass in the definition of kinetic en- 
ergy). Here F{9,p) is the normalized distribution func- 
tion of the system. Eq. ([7]) furnishes the link between 
Tjfin and the thermodynamical temperature T^h = 
For a Boltzmannian system: F{9,p) = Cexp(— /3p^/2), 
where C is a normalization constant [l| . Substituting the 
latter expression into Eq. ^ one gets Tkm = !//? = Tth- 
In this case the statistical and kinetic definitions of spe- 
cific heat provide the same result: (ct,)^j^ = {cv)i^i^- 

From the above, one can conclude that, for a Boltz- 
mannian system, the canonical specific heat is always 
positive, regardless the definition of temperature. As a 
consequence, for a Boltzmannian LRI system, the pres- 
ence of negative specific heat in the microcanonical en- 
semble is the signature of ensemble inequivalence. 

A completely different scenario may arise in the case 
of non-Boltzmannian or out-of-equilibrium LRI systems. 
As we shall demonstrate, for those systems one can ex- 
pect to measure negative specific heat in the canonical 
ensemble. 

Consider the "fermionic" system characterized by the 
entropy paradigmatic example of a LRI system 

"trapped" in a QSS. Following the procedure outlined 
above, one can perform the LET of s[f] and calculate 
the rescaled free energy Requiring the free en- 

ergy to be stationary, one gets the following distribution: 



fQSs{0,p) = 



fo 



1 + ef^fo^P - ^MfQSs] cose- A/,[/qss] sinfl) + a 

(8) 

Here a plays the role of a Lagrange multiplier associ- 
ated, rispectively, to mass conservation, while M^lf] = 
J fcos{0)d0dp, My[f] = J fsm{0)d0dp stand for the two 
components of magnetization in the N oo limit. It is 
worth stressing that the expression ([H]) for the stationary 
distribution in the canonical ensemble is formally identi- 
cal to the one obtained by requiring the entropy (O to 
be stationary in the framework of the microcanonical en- 
semble 0. In that case, /3 and a are Lagrange multipliers 
associated to energy and mass conservation. 

Fig. [1] shows the entropic curve, s = s(e), for a given 
initial condition /q. The curve is calculated in the frame- 
work of the microcanonical ensemble, using the following 
procedure. Once fixed e and /o, one first determines the 
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FIG. 1: Lynden Bell entropic curve, s = s(e), calculated for 
the initial condition fo = 1.1096, in the framework of the mi- 
crocanonical ensemble. Inset: Difference, As, between mag- 
netized and non-magnetized branches of the entropic curve, 
as a function of energy. A chan ge o f sign in As(e) is the signal 
of a phase transition, see also [l5|]. There is a second-order 
phase transition aX e — 0.6123 and a first-order phase transi- 
tion at e = 0.5987. Note that the latter cannot be noticed on 
the (too large) scale used to plot the entropic curve. 



corresponding values of , My , a and /? in the expres- 
sion ([8|) for /qss- This can be done by imposing energy 
and mass conservation and by using the definition of 
and My Q. In general, one obtains several possible dis- 
tributions /qss, corresponding either to magnetized (i.e., 

M = {M^ + M^y^^ ^ 0) or to non-magnetized (i.e., 
M ~ 0) states. The value of s(e) is found by substituting 
the possible solutions into ^ and by solving numerically 
the integral: the retained value is the one corresponding 
to maximum. This also selects the expected (magnetized 
or non-magnetized) stable equilibrium distribution. 

The Lynden-Bell entropic curve s(e) is always concave 
(see Fig. [T]): such a result is a clear indication of ensem- 
ble equivalence. Assuming F{0,p) = /qss(^:P) in Eq. 
(O, one gets a formal relation between the kinetic tem- 
perature, Tkin, and the Lynden-BcU Lagrange multiplier 
/? = 1/Tkin (inverse thermodynamic temperature). 

Reconsidering now equations ([5]) and with T = 
Tkin, one is led to the surprising conclusion that the 
measured specific heat can be negative both in the mi- 
crocanonical and in the canonical ensembles, no matter 
the convexity of the entropy and of the free energy func- 
tions. In fact, in both cases the sign of the specific heat 
also depends on that of the derivative of the function 

The curve Tkin = TkiniP) can be calculated in the 
framework of the canonical ensemble, as follows. Once 
fixed (3 and /o, one determines the stable equilibrium dis- 
tribution by finding out the corresponding values of Mx, 
My and a. This can be done by imposing the normaliza- 
tion constraints and by using the definition of magneti- 
zation The obtained result is then substituted into 
Eq. ([7|) and the integral solved numerically. The curve 
Tkin = 7kin(/3) is shown in Fig. [Da,. 
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FIG. 2: Fig. a): Kinetic temperature, Tkin, as a function 
of the Lynden-Bell Lagrange multiplier /3, calculated for the 
initial condition fo = 1.1096. Fig. b): Kinetic caloric curve, 
l/7kin vs. e (continuous line), and thermodynamic caloric 
curve, 1/Tth — P vs. e (dotted line), calculated at /o — 
1.1096. 

As it can be seen, one gets dTkin/idp) > and, thus 
negative specific heat, {cv)^^^^ < 0, in both ensembles, in 
the magnetized region. The change of sign occurs in cor- 
respondence of phase transitions: from positive to nega- 
tive in correspondence of a second-order phase transition, 
at smaller (3 values, and from negative to positive in cor- 
respondence of a first-order phase transition, at larger /3 
values, where the curve displays a discontinuity. Note 
that this behaviour is fundamentally different from that 
displayed by Boltzmannian LRI systems. Indeed, as al- 
ready mentioned, in that case the change from positive 
to negative (thermodynamical) specific heat may occur 
only in correspondence of a first-order phase transition 
[W^ . Phase transitions are also signaled by a change of 
stability of the magnetized and non-magnetized branches 
of the entropic curve, as shown in the inset of Fig. [TJ The 
presence of a negative kinetic specific heat in the mag- 
netized region can be also noticed by looking at the ki- 
netic caloric curve shown in Fig. [5]d (continuous line). 
As expected on the basis of the result shown in Fig. 
m the thermodynamic caloric curve (dotted line in Fig. 
[2}j) , is instead characterized by a negative derivative over 
the whole energy range, corresponding to (c„)j.[^ > in 
both ensembles. It is worth stressing that the occurrence 



of negative specific heat in the microcanonical ensemble 
when the system is magnetized is confirmed by the results 
of direct A^-body simulations, based on the Hamiltonian 

As we anticipated, a further surprising result is that 
the presence of negative specific heat is associated to 
equivalent statistical ensembles. A clear indication of 
ensemble equivalence can be obtained by inspection of 
Fig. [U showing a concave entropic curve, and of Fig. [2Ji, 
showing the occurrence of positive and negative canonical 
and microcanonical specific heats in the same interval of 
/? values. The same indication can be drawn by the ther- 
modynamical caloric curve shown in Fig. ^jp. However, 
the curves displayed in figs. [T] and [2] have been obtained 
for a particular initial condition /q. In order to investi- 
gate ensemble equivalence for any fo, we have compared 
the values of M obtained by solving the Lynden-Bell 
problem in the two ensembles for different initial con- 
ditions, following the above outlined procedures. Results 
(not reported here) show that the magnetized and un- 
magnctized regions perfectly coincide, and this confirms 
ensemble equivalence. 

Our study demonstrates the unexpected possibility to 
measure negative specific heat in the statistical canoni- 
cal ensemble. As we argued above, this conclusion docs 
not contradict any fundamental law of thermodynamics 
but ultimately originates from the non-Boltzmannian fea- 
tures that are possibly associated to non-equilibrium, as 
well as equilibrium, dynamics. With reference to the case 
at hand, and for a specific window of the parameter space 
that we traced back to magnetized QSS, the particles 
gain in kinetic energy when an energy quota is passed 
from the system to the heat reservoir. This surprising ef- 
fect arises spontaneously and it is driven by the inherent 
ability of the system to self-organize at the microscopic 
level. The average particle velocity is enhanced at the 
detriment of the potential energy, and this yields in turn 
to an increase of the experimentally measurable kinetic 
temperature. The opposite holds if the energy flows to- 
wards the system: the particles cool down, while the po- 
tential contribution grows so to guarantee for the needed 
energy balance. This phenomenon, that we have here 
demonstrated with reference to the HMF model, is in 
principle general and can potentially extend to all those 
settings where non-Boltzmannian effects play a role. It 
is particularly attractive to speculate on the possibility 
of realizing efficient thermal devices, working with a non- 
Boltzmannian fluid. Can one improve over current im- 
plementation by exploiting the aforementioned tendency 
to concentrate the residual energy amount into a kinetic 
quota that supports and enhances the particle motion? 
As a final comment, let us emphasize that, by invok- 
ing the above mentioned analogy with the HMF model, 
dedicated CARL experiments should allow for a direct 
verification of the predictions here elaborated and so con- 
tribute to shed light on the proposed scenario. 
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